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Systematic variation of the stellar Initial Mass Function with 
velocity dispersion in early-type galaxies 



(N 

o 

(N 
^— > 

O 

O 



o 

u 

6 

(N 
> 

in 
o 

(N 



13 



Ignacio Ferreras 1 , Francesco La Barbera 2 , Ignacio G. de la Rosa 3 ' 4 , Alexandre Vazdekis 3 4 , 
Reinaldo R. de Carvalho 5 , Jesus Falcon-Barroso 3 ' 4 , Elena Ricciardelli 6 

1 Milliard Space Science Laboratory, University College London, Holmbury St Mary, Dorking, Surrey RH5 6NT 

2 INAF - Osservatorio Astronomico di Capodimonte, Napoli, Italy 

3 Institute de Astrofi'sica de Canarias, C/ Via Ldctea s/n, La Laguna, E-38200 La Laguna, Tenerife, Spain 

4 Departamento de Astrofi'sica, Universidad de La Laguna, E-38205 La Laguna, Tenerife, Spain 

5 Instituto Nacional de Pesquisas Espaciais/MCT, S.J. dos Campos, Brazil 

6 Departament d'Astronomia i Astrofisica, Universitat de Valencia, C/Dr Moliner 50, E-46100, Burjassot, Valencia, Spain 



MNRAS Letters: Accepted 2012 October 11. Received 2012 October 11; in original form 2012 August 15 



ABSTRACT 

An essential component of galaxy formation theory is the stellar initial mass function (IMF), 
that describes the parent distribution of stellar mass in star forming regions. We present ob- 
servational evidence in a sample of early-type galaxies (ETGs) of a tight correlation between 
central velocity dispersion and the strength of several absorption features sensitive to the pres- 
ence of low-mass stars. Our sample comprises ~ 40,000 ETGs from the SPIDER survey 
(z <, 0.1). The data - extracted from the Sloan Digital Sky Survey - are combined, reject- 
ing both noisy data, and spectra with contamination from telluric lines, resulting in a set of 
18 stacked spectra at high signal-to-noise ratio (S/N }t 400A _1 ). A combined analysis of 
IMF-sensitive line strengths and spectral fitting is performed with the latest state-of the art 
population synthesis models (an extended version of the MILES models). A significant trend 
is found between IMF slope and velocity dispersion, towards an excess of low-mass stars in 
the most massive galaxies. Although we emphasize that accurate values of the IMF slope will 
require a detailed analysis of chemical composition (such as [a/Fe] or even individual element 
abundance ratios), the observed trends suggest that low-mass ETGs are better fit by a Kroupa- 
like IMF, whereas massive galaxies require bottom-heavy IMFs, exceeding the Salpeter slope 
at cr > 200kms" 1 . 

Key words: galaxies: elliptical and lenticular, cD - galaxies: evolution - galaxies: formation 
- galaxies: stellar content - galaxies: fundamental parameters. 



1 INTRODUCTION 

The Initial Mass Function (hereafter IMF) is denned as the distri- 
bution of stellar masses in a single population at the time of birth. It 
has been usually considered a universal function, partly because of 
the complexities in obtaining proper observational constraints. The 
traditional approximation by a single power law ( Salpeter] EUl) 
has undergone numerous updates, with more complex functions 
that include a significant flattening of the slope towards low-mass 
stars dMiller & Scalo Ill979l; IScalo II 19861; iKroupa Il200ll ; IChabrierl 
2003). For a recent revie w on the IMF and its possible variations, 
see Bastian et al. 

Past analyses based on dynamical modelli ng favoured 
Kroupa-like IMFs in ETGs JCappellari et al. 1 120061) . This result 
was also found wit h strong gravitation al lensing on the bulge of 
a late-type galaxy JFerreras et al. 1 1201 01) . However, the IMF was 
suspected to depend on galaxy mass, when extending the analy- 
sis to more massive systems, as in the dynamical modelling of 



Thomas et al. with galaxy morphology also playing a 

role jBrewer et al. 1120121) . Constraining the low-mass portion of 
the IMF is a challenging task, as faint, low-mass stars do not 
contribute significantly to the integrated spectrum. The first ob- 
servational atte mpts to constr ain th e low-mass end of the IMF 
were made by ICohen I dl978h and iFaber & French] dl980b , to- 
wards the centers of M31 and M32, measuring several spectral 
features sensitiv e to the rela t ive co ntribution of M dwarves and 
M giants. Later, ICarter et"al! (1 1986b extended the study to a sam- 
ple of massive early-type galaxies, and found that Nal was en- 
hanced, especially in massive galaxies, with strong radial gradi- 
ents (being most enh anced in the central region). More recently, 
Cenar ro et al. 1 d2003l) proposed a relation towards an excess of low- 
mass stars in massi ve galaxies, from a study of the Call triplet 
region at ~8500 A . IVan Dokkum & Conroy] j2010l) probed two 
features in the NIR that are strongly dependent on the fraction of 
low-mass stars, to find a bottom-heavy IMF in a sample of four 
Virgo ETGs. This result was recently extended to an additional set 
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Figure 1. Trend of the three absorption features targeted in this paper (TiO 1 , 
Ti02 and Na8190) with respect to central velocity dispersion, for stacked 
SDSS spectra of ETGs from the SPIDER sample. The SEDs have been 
smoothed to a common velocity dispersion of 300 km s — 1 with the spectral 
resolution of SDSS, and continuum subtracted with a second order polyno- 
mial. We only show, for clarity, three out of the 1 8 stacks, spanning the full 
range of velocity dispersion. The shaded regions mark the positions of the 
red and blue sidebands, and the central bandpass. 
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Figure 2. Trend of the equivalent widths of TiOl, Ti02 and Na8190, with 
respect to velocity dispersion. All measurements are performed on data con- 
volved to a common velocity dispersion of 300 kms -1 with the spectral 
response of SDSS. The error bars give the statistical error from the stacks, 
at the 3 <r level. The shaded regions on the right correspond to model pre- 
dictions for three choices of IMF unimodal slope, as labelled, spanning a 
range of ages from 5 to 10 Gyr and metallicity from log Z/Zq = —0.4 to 
+0.2 



of 34 ETGs from the SAURON survey, with similar conclusions 
(Van Dokkum & Conrov"ll2012h . favouring the interpretation of a 
non-universal IMF. An alternative scenario to explain these data 
wo uld invoke over/ under abundances of individual elements (see 
e.g. lWorthevlll998h . Bottom-heavy IMFs would imply that mas- 
sive galaxies observed at high redshift are even more massive than 
the standard assumption of a Kroupa-like IMF, posing challenging 
constraints on the physics underlying the conversion of gas into 
stars during the first stages of galaxy formation. This issue is of 
special relevance to the standard paradigm of hierarchical galaxy 
formation, whereby massive galaxies are supposed to grow mainly 
through the merging of sma ller systems, with the growth rate lim- 
ited by cosmology (see, e.g jDe Lucia et al. l2006h , requiring alter- 
native chan nels of massive galaxy f ormation, such as cold accretion 
(see, e.g. lBirnboim & Dekel 1120031) . 

Simple theories of star formation suggest a power law 
behavio ur, with a truncation at the low-ma ss end, around 
0.3 Mn dLarsonll2005l : [McKee & Ostriker 2007), consistent with 
some of the most popular IMFs used dKroupal l200ll ; IChabrierl 
l2003h . However, due to the lack of strong observable tracers in the 
observed spectra arising from low mass stars, a robust constraint at 
the low-mass end of the IMF has remained elusive for years. In ad- 
dition to the work cited above, more recent evidence has been ac- 
cumulating dVan Dokkum & Conrov 1120101 : ISpiniello et al. 112012k 
ISmith et al. Il2012f) . towards the possibility of a non-universal IMF. 
This hypothesis is further supported by an independent study that 
made us e of detailed dynamica l modelling of nearby spheroidal 
galaxies I Cappell ari et al. Il2012l) . Strong gravitational lensing over 
galaxy scales also hint at bottom heavy IMFs in massive galaxies 
dAuger et al. 1120101) . The implications of a change in the IMF war- 
rant further analysis of this problem. 



2 SAMPLE 

The hypothesis of a systematic change in the Initial Stellar Mass 
Function with respect to global properties of galaxies, such as 
mass or luminosity, are optimally tested on samples of galaxies 
that harbour stellar populations as homogenous as possible, for 
a clean analysis of the spectra. The SPIDER sample (Spheroids 
Panchromatic Investiga tion in Different Environmental Regions, 
lLa Barbera et al. l|2010l) is optimal for this purpose, as it comprises 
39,993 nearby (0.05<z<0.095) ETGs, in a wide range of masses , 
from 2 x 10 10 to 10 12 Mrr, (dynamical mass. lde la Rosa et al. |2012|) . 
Spectra of these galaxies were ret rieved from the sixth da ta release 
of the Sloan Digital Sky Survey dAbazaiian et~al~ll2009h (SDSS), 
de-redshifted to a common rest-frame (with 1 A binning) and cor- 
rected for foreground Galactic extinction (see lLa Barbera et al. I 
[201Cj, for details). In order to test variations of the IMF, we as- 
semble 18 stacked spectra spanning over a wide range of central 
velocity dispersions (a) between 100 and 320 kms -1 in steps of 
10km s^ 1 , except for the two bins for the highest values of a, 
where - because of the smaller number of galaxies in the sample 
- we adopted the range [260,280], and [280,320] km s _1 , respec- 
tively. We select the sample according to a, because (i) the under- 
lying stellar populatio ns of ETGs are know n to correlate strongly 
with this observable teemardiet"ai~ll2005l) . and (ii) the spectral 
analysis depends sensitively on velocity dispersion, as a higher a 
introduces an effectively lower spectral resolution, an issue that we 
take into account for a consistent comparison (see below). 

Given the redshift range of the SPIDER sample, the 3 arcsec 
diameter of the SDSS fibers maps the central ~ 2 — 3.5 h^ 1 kpc. 
For each bin, we perform median stacking of the available spectra, 
considering only pixels with no SDSS flag raised. In order to avoid 
possible biases related to intrinsic dust and differences in signal to 
noise ratio (S/N) within the bin, we excluded spectra in the lowest 
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quartile of the S/N distribution, as well as those galaxies (14% of 
the entire sample) with E(B — V) > 0.1 mag. The colour excess 
is measured by fitting each SPS S spectrum with a variety of stellar 
population models, assuming a Cardelli, Clayton. Mathis 1 Jl989l) 
extinction law, as detailed in Swin dle et alj 1 201 1 ). 



Spectral Fit Na81 



3 IMF-SENSITIVE INDICATORS 

The SDSS spectral range allows us to study several IMF-sensitive 
spectral features, such as the Nal doublet at AA8183,8195A 
jSchiavonet"aLlll997h . This feature is prominent in the atmo- 
spheres of low-mass dwarf stars. Therefore, its strength in a stellar 
population depends quite sensitively on the shape of the IMF at the 
low-mass end. Over the IMF-sensitive wavelength range (8183- 
8 195 A), our stacked spectra feature a remarkably high S/N, with 
an average of ~ 400 per A for the lowest and highest a bins, 
reaching values as high as S/N ~ 850 per A in the stack with 
velocity dispersion a ~ 160 kms -1 . We use here a slightl y mod- 
ified v ersion of the NaI8200A index proposed by Vazdeki s et al. I 
< l2012h . defined as Na8190, where the blue passband is measured 
in the [8137,8147]A interval (instead of [8164,8173]A), where the 
MIUSCAT models are found to give a better match to the data, 
within ~ 0.3% (compared to 1% at A ~ 8170A). The red and 
central bandpasses are the same as for the NaI8200A index. 

We performed extensive tests to check for possible contam- 
ination of the Na8190 feature from sky lines, by combining (1) 
all galaxies at each A regardless of the SDSS flag value; (2) only 
flux values not affected by sky absorption (i.e. lA apart from any 
telluric line) and sky emission; and (3) only ETGs at the lowest 
redshifts in the SPIDER sample (0.05 < z < 0.07), where the Nal 
doublet is observed in a region free from any significant sky con- 
tamination. For all a bins, the variation of the average Nal flux, 
using the above three stacking procedures, was found to be smaller 
than 0.06%, i.e. negligible for our purposes. 

In order to increase the robustness of the analysis against pos- 
sible biases regarding an overabundance of [Na/Fe] in massive 
galaxies, we also consider two line strengths targeting TiO, a fur- 
ther discriminant with respect to the presence of low-mass stars. 
We use the TiOl and TiQ2 defi nitions from the standard Lick sys- 
tem (see e.g. Trager et al. 1998). The Ti02 index has been recently 
used by ISpiniello et al. I d2012r) to constrain the IMF of a num- 
ber of early-type gravitational lenses. Fig. [Tj shows three of our 
18 stacks, in the regions of interest for this work. All SEDs have 
been smoothed to a common velocity dispersion of 300 km s ~~ 1 for 
a consistent comparison, and were continuum subtracted, using a 
second order polynomial. A significant trend is apparent from low- 
to high- velocity dispersion of all indices. This trend encodes a com- 
plex range of variations in the underlying stellar populations, most 
notably age, metallicity, abundance ratios, and IMF. 

Note in the Na8190 region (bottom), the centroid of the line 
shifts bluewards as a increases. However, because of the asymme- 
try of the spectral region around the line, more flux would be ex- 
pected to contribute from the red side of the line, therefore shifting 
the position of the centroid redwards. This shifting reveals that it is 
the absorption of Na at 8 190 A - and not other contaminating lines 
in the vicinity - that changes with respect to the velocity dispersion 
of the ETG. The two most obvious interpret ations of this eff ect in- 
volve either an overabundance of [Na/Fe] dWorthevlll9"98h , or a 
change in the IMF dVan Dokkum & Conrov 1120101) . In a forthcom- 
ing paper (La Barbera et al., in preparation) we explore in detail 
the effect of individual overabundances. However, in this letter, we 
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Figure 3. The 68, 95 and 99% confidence levels of the likelihood from 
individual constraints (from left to right) of spectral fitting, Na8190, TiOl 
and Ti02 are shown on the age-metallicity plane for two of the stacks, as 
labelled by their central velocity dispersion. Bimodal IMFs are used here, 
although the results for the unimodal case are similar. 



combine the three line strengths - which rely on different species - 
to confirm the trend towards a bottom-heavy IMF in massive galax- 
ies. 

Fig. E] shows the line strengths of the stacked SEDs, for 
the three IMF-sensitive indices. The strengths are corrected to 
a common broadening of cr re f =300kms~ 1 (in addition to the 
wavelength-dependent SDSS resolution). A strong correlation of 
these line strengths with velocity dispersion is evident. We checked 
that the fixed aperture of the fibres used by SDSS to retrieve spec- 
tra did not introduce a bias with respect to size, by comparing the 
change in line strengths with respect to the angular extent of the 
effective radius. The variation of the indices, at fixed a, is around 
an order of magnitude smaller than the trend shown in Fig. [2] The 
shaded regions on the right of the figure motivate the methodol- 
ogy followed in this paper. Each one gives the model predictions 
for a choice of unimodal IMF slope (labelled), over a range of SSP 
ages and metallicities. The values of all three indices at high ve- 
locity dispersion can only be reconciled with a bottom-heavy IMF. 
Our methodology - explained below - consists of removing the 
degeneracies from age and metallicity by combining line strength 
information with spectral fitting over the optical range. 



4 CONSTRAINING THE IMF 



We make use of MIUSCAT IVazdekis et all |2012| ; 
iRicciardelli et al. I l2012h . the spectrally extend ed version of 
the s tellar population synthesis models MILES (Vazdek is et al. I 
2010), in order to map the systematic trend of the TiOl, Ti02 and 
Na8190 line strengths with respect to IMF slope. These models 
combine state-of-the art stellar libraries in the optical and NIR 
windows, creating a set of spectra with a uniform resolution of 
2.51 A throughout the wavelength range 3465-9469 A . Our data 
are compared with grids of simple stellar populations (SSPs), 
assuming either a u nimodal or a bimodal IMF (as defined in 
IVazdekiset"al~i ri996). over a stellar mass range 0.1 — 100 Mq, 
For the unimodal case, we adopt the logarithmic slope V = x — 1, 
wh ere dN/dm oc m~ x is the Initial Mass Function, such that 
the ISalpeterl i ll 9551) IMF corresponds to Y = 1.35. A bimodal 
IMF replaces the M < 0.6Mq interval by a flat portion and 
a spline to match the power law at the high mass end. It gives 
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a closer representati on of Kroupa-like IM Fs for F = 1.3 (see, 
e.g. Figure 1 of IVazdekis et ail 120031) . We use MIUSCAT 
SSPs spanning a wide range of ages (1-13 Gyr) and metal- 
licities (-1.0 < log Z/Z e +0.22), for different values of 
T = {0.3,0.8, 1.0, 1.3, 1.5, 1.8,2.0,2.3,2.8,3.3}. In order to 
obtain an accurate estimate of a for each sta ck, we perform spec- 
tral fi tting with the software STARLIGHT dCid Fernandes et al. I 
120051) . in the range 3800-8400 A. STARLIGH T can be used to 
extrac t full star formation histories (see e.g. Ide la Rosa et al. I 
12012b . but we are only interested here in assessing the robustness 
of the measured a with respect to the basis SSPs. We find no 
significant trend when choosing template SSPs with different 
IMFs, with a variation Act ^ 1 kms -1 . Furthermore, for each 
stack, the a determined by STARLIGHT is consistent, within 
~ 10%, with the median a of the stacked spectra from SDSS. 

We parameterise a stellar population by an age, metallic- 
ity, and IMF slope. More detailed models will be presented 
elsewhere (L a Barbera et a l., in preparation). We note that, 
according to iThomas et al, ] feOllah . both TiOl and Ti02 arc 
rather insensitive to non-solar abundance ratios. Furthermore, 
IConrov & van Do kkum find an anticorrelation between 

[a/Fe] and NaI0.82 strength (similar to our Na8190 index). Hence, 
we do not expect our results, regarding the non-universality of the 
IMF, to be affected by the well-known overabu ndance of [a/Fe] 
in massive galaxies (see e.g. iTrager et"aT, 1 l2000l) . To confirm this 
point, we studied two subsamples within the highest-cr bin, split 
with respect to the distribution of [a/Fe] values, constructing 
low-(~ +0.1) and high-(~ +0.3) [a/Fe] stacks. The measured 
Ti01(mag)/Ti02(mag)/Na8 190(A) strengths are 0.0351 ± 0.0014 
/ 0.0860 ± 0.0010 / 0.74 ± 0.07, for the high [a/Fe] subsample, 
0.0352 ± 0.0015 / 0.0833 ± 0.0010 / 0.77 ± 0.08, for the low 
[a/Fe] subsample, which is fully consistent with the original stack, 
including all values of [a/Fe]: 0.0364 ± 0.0006 / 0.0842 ± 0.0005 
/ 0.78 ±0.03. 

Our analysis follows a hybrid method, whereby two indepen- 
dent probability distribution functions for age, metallicity and IMF 
slope are obtained using spectral fitting and line strengths, respec- 
tively. For the former, we fit the region 3900-5400A , using SSPs 
from MIUSCAT, after convolving them from the orig inal 2.51 A 
resolution, to the fiducial velocity dispersion of 300 kms -1 , plus 
SDSS spectral resolution. A grid of 32 x 32 x 10 models in age, log 
Z and r is used to obtain a probability distribution function (PDF), 
via the likelihood C oc exp(— Ax 2 /2). For each line strength, we 
also define a second PDF in the same grid. The joint PDF is defined 
as the product of the corresponding (independent) PDFs. Finally, 
we marginalize over all parameters but F, to obtain the PDF of the 
IMF slope. Fig.[3]illustrates the constraining power of the different 
observables, individually, on the age and metallicity of the model 
SSP. Notice that the combined analysis of the hybrid method is ac- 
ceptable for all three cases (i.e. the joint likelihood does not reduce 
to the product of mutually unlikely regions of parameter space). 



5 CONCLUSIONS 

Fig. [2] motivates the need to invoke a non-universal IMF when 
comparing galaxies over a range of velocity dispersions, with 
an increased contribution from low-mass stars in the most mas- 
sive galaxies. Fig. [4] presents the best-fit slope of the IMF, when 
marginalizing over age and metallicity, as a function of velocity dis- 
persion, given by the probability weighted estimates using the joint 
PDF, for three cases, depending on whether TiOl (filled circles), 
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Figure 4. Variation of the IMF slope - unimodal (top) and bimodal (bot- 
tom) distributions - against central velocity dispersion. The shaded region 
corresponds to the 68% confidence level of the joint PDF including spectral 
fitting and all three line strengths (TiOl, Ti02 and Na8190). The Salpeter 
(unimodal) and Kroupa-equivalent (bimodal) cases are shown as horizon- 
tal dashed lines. A least squares fit to the data is shown for reference (solid 
line). The rightmost panels give the stellar mass-to-light ratios in the SDSS- 
r band for a 10 Gyr old population at solar metalliticy, illustrating the large 
variations one could expect depending on the choice of IMF. 

Ti02 (open squares), or Na8190 (crosses) are used in the analysis. 
Both unimodal (top) and bimodal (bottom) IMFs are considered. 
We note that at low velocity dispersion, the constraint on IMF slope 
becomes more challenging, as l ow-mass ETGs feature complex 
star formation histories (see e.g. Ide la Rosa et al. I [20121) . making 
the SSP approach used here not fully applicable. Therefore, Fig. [4] 
only shows the range 150-300 km s - , where the approximation of 
a SSP is justified. The general trend towards a bottom-heavy IMF 
is evident in all cases. The result for the joint PDF corresponding to 
all three indices plus spectral fitting is shown as a grey shaded re- 
gion - extending over the 68% confidence level - with a black line 
giving a simple least squares fit to the data. We emphasize here that 
accurate values of the IMF slope will require a detailed analysis of 
abundance ratios (such as, e.g., [a/Fe] or [Na/Fe]). The predicted 
IMF slopes for massive ETGs do not pose a problem to optical-NIR 
photometry, where the contribution from low-mass stars would be 
most important. As an example, the MIUSCAT models for a 10 Gyr 
population at solar metallicity give V — K = 3.04 for a (bimodal) 
IMF slope T = 0.30, versus V - K = 3.03 for T = 2.80 (A 
unimodal distribution will give an excess in V — K of ~ 0.25 mag 
over a similar range of F). Hence, a variation in the IMF is com- 
patible with the observed V — K colours. Furthermore, V — K 
colours alone cannot be used to constrain the IMF, a well-known 
result, that explains why constraints on the IMF have remained 
elusive. The V — K colour of our sample ranges from 3.0 at the 
lowest velo city dispersion bin, to 3 .3 at a velocity dispersion of 
300 km s^ dLa Barbera et al. lioich . 

The rightmost panels of Fig. [4] should serve as a warning 
to applications of these trends to infer stellar masses. The M/L 
in the SDSS-r band is shown for a typical old, metal-rich popu- 
lation, such as those found in these galaxies. The difference be- 
tween a unimodal (top) and a bimodal (bottom) IMF is very sig- 
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nificant for massive galaxies, whereas the fits for either choice of 
IMF are equally acceptable. The correlation presented in this let- 
ter is consistent with studies based on different methods, involv- 
ing strong gravita t ional lensing, or dynamical modelling (see, e.g . 
Ifreu et al. I lioiol; iThomas et all 1201 lbt ICappellari et all 120121) . 
Regarding the st udy of stellar populati ons, our work is consistent 
with the trend in lCenarro et al. I d2003l) . based on a spectroscopic 
analysis of the Ca ll triplet region. In addit ion, the recent photo- 
metric analysis of iRicciardelli et al71 (2012, see their Figure 16), 
is compatible with our a nalysis. Our results are a lso in agree- 
ment with the findings of IVan Dokkum & Conrovl d2010l l2012h . 
We emphasize that our analysis is based on a completely inde- 
pendent stellar library to their work. For instance, in the region 
around the Na8190 feature, the MIUS CAT models u se stellar spec- 
tra from t he Indo-US library dValdes et al. 1 12004b . whereas the 
models of Conrov & van Dokkum (2012) use the IRTF library 
dCushing et al. ll2005l ; lRavner et al. 12 009 ) . In a forthcoming paper, 
we will explore in detail different aspects related to the method- 
ology, including the contribution from enhacements of individual 
elements, non-solar [a/Fe], or composite stellar populations. The 
large size of our sample - comprising nearly 40,000 ETGs - al- 
lowed us to go beyond a simple test of non-universality of the IMF, 
enabling us to obtain a trend between the IMF slope and veloc- 
ity dispersion in more detail than previously found. The existing 
correlation between V and velocity dispersion in ETGs suggests a 
significant transition in the properties of star forming regions as a 
function of galaxy mass, an issue with important implications on 
galaxy formation theories. 
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